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1. Introduction 



In continuum quantum field theory in four dimensions, massless Dirac fermions respect chiral 
symmetry. This symmetry is crucial for our understanding of low energy processes of QCD: since 
the light quark masses are small compared to the QCD scale, one would expect to see signatures 
of chiral symmetry in the spectrum. Since these signatures are not observed, we assume that QCD 
dynamically breaks chiral symmetry at some scale A^; this leads to concrete predictions which 
are phenomenologically successful. Consequently, one would like latticized QCD to enjoy some 
analogue of chiral symmetry. This is achieved if the quarks satisfy the Ginsparg-Wilson relation [|l|] ; 
current realizations are domain wall fermions and overlap fermions. These quark discretizations 
are collectively known as Ginsparg-Wilson (GW) quarks. 

Unfortunately GW quarks come with a price. Simulations using GW quarks are slower [§] than 
simulations with Wilson or staggered quarks. Staggered quarks, in particular, have been popular 
recently due to the availability of the MILC lattices [^]. But this discretization scheme has its own 
issues: staggered chiral perturbation theory [|], ^ |6|] is theoretically complicated, as demonstrated, 
for example, by the large number of operators in the NLO Lagrangian [f^]. 1 

A practical compromise between the GW and staggered discretization schemes is to use stag- 
gered fermions for the sea quarks and GW fermions for the valence quarks [||]. This is particularly 
attractive since one can employ the MILC lattices, but as we shall see, it is also theoretically clean. 
Several recent simulations [||, [K], 11, 12, 13, 14] have utilized this kind of mixed action scheme, 



so it is important to understand the chiral extrapolation appropriate to these simulations. 

Nevertheless, these mixed action simulations face their own challenges. Chiral perturbation 
theory for mixed action simulations (MA^PT) [15, [If]] depends on a new, as yet undetermined 



parameter, Cmi x - This parameter shifts the mass of mixed valence-sea mesons at tree level, so that 
for a meson with valence quark v and sea quark s, 

m 2 vs = B(m v + m s )+ a 2 A M ix, (1.1) 

where Ajviix = 16CMix// 2 and a is the lattice spacing. In addition, a mixed action theory violates 
unitarity at any finite lattice spacing, so one must take the continuum limit before unitarity can be 
restored, unlike a partially quenched theory. We will comment on these issues below. 

Before we continue, let us introduce some notation. We shall denote our valence quarks as u,d 
and s while our sea quarks are j,l and r. The leading order (LO) mass of a meson of quark content 
a,b will be denoted m a [,. Meson masses which suffer a tree level lattice correction will be denoted 
with a tilde; for example, the mass of a meson with quark content j,j is mjj = 2Bmj + c? 2 A/, where 
a 2 Aj is the mass shift of a taste identity meson due to the staggered potential. The next-to-leading 
order (NLO) pion mass will be denoted m n and similarly for the kaon mass. The leading order 
decay constant is denoted by / while the NLO pion and kaon decay constants are f n and fg, 
respectively. 



1 We assume the validity of the fourth-root trick, though it has not been shown that the continuum limit of latticized 
QCD with rooted-staggered quarks is in the same universality class as QCD. 
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2. NLO Corrections 



Recently, the %% scattering length was computed [ |17| ] in mixed action chiral perturbation 
theory (MA^PT) for GW valence quarks on a staggered sea. In terms of leading order parameters, 
the scattering length is given by 
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mjj — m n . This result simplifies considerably [ ]T7| ] if one rewrites the scattering length 



where Apg 

in terms of the mass of the pion on the lattice and the decay constant of the pion measured on the 
lattice (we refer to these as lattice-physical parameters.) Then the scattering length becomes 
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It is worth commenting on some aspects of this simplification. In Eq. (|2.2|), the unknown coun- 
terterm which must be determined on the lattice, l n]t , is identical to the continuum, unquenched 



counterterm at NLO, while in Eq. (2.1) there are additional unknown unphysical terms 1' P q and l' ul . 
These unphysical effects must somehow be removed in the analysis of the lattice data. The expres- 
sion Eq. ( |2.1[ ) depends on the masses of the mixed valence -sea meson, nij u . This mass depends 
on the unknown parameter Cmi x ; therefore, this parameter must be determined if one wishes to 



extrapolate using Eq. ( |2.1| ). The lattice-physical expression Eq. ( |2.2| ), however, does not depend on 
this unknown parameter; in fact, the lattice-physical expression does not depend in any way on the 
mixed valence-sea mesons, and only differs from the continuum, unquenched scattering length by 
a known term which has its origin in the unitarity violating flavour neutral propagators. Our goal in 
this section is to explain why these simplifications occur, and to what extent we may expect similar 
simplifications in other mesonic processes. 

To begin, let us consider the structure of the NLO Lagrangian of MA^PT. This Lagrangian 
is determined by the symmetry structure of the underlying partially quenched, mixed action form 
of QCD, supplemented by the assumption that chiral symmetry is spontaneously broken. Since 
the theory has two sectors — the valence and sea sectors — it is worth considering the spurions 
arising from these sectors separately. In the valence sector, the quarks satisfy the Ginsparg-Wilson 
relation; therefore, they only explicitly violate chiral symmetry through the quark mass term. This 
is the same as in continuum, unquenched QCD and therefore the purely valence sector of the 
NLO chiral Lagrangian in the mixed action theory is the Gasser-Leutwyler Lagrangian. There 
are additional sources of chiral symmetry violation in the sea sector. In particular, for staggered 
sea quarks, there are additional spurions associated with taste symmetry violation at finite a. One 
might expect that these spurions would lead to a dependent analytic terms in the NLO scattering 
amplitudes; however, we will now show that this does not occur. 

At NLO in the power counting scheme appropriate to current simulations, the NLO operators 
in the chiral Lagrangian only contribute at tree level. Of course, all of the in and out states in a 
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simulation involve valence quarks, so for the purposes of an NLO computation one can set the sea 
quarks equal to zero in the NLO operators. Since the spurions arising from the sea sector go to zero 
when the sea quarks go to zero, we see that at NLO the only relevant spurions are the purely valence 
spurions. However, we have seen above that the NLO purely valence Lagrangian coincides with 
the Gasser-Leutwyler Lagrangian, so no lattice spacing dependent analytic terms can arise in this 
way. An exception to this argument occurs for double-trace operators; however, these operators can 
only renormalize the parameters / and B occurring in the LO chiral Lagrangian. In lattice-physical 
parameters we eliminate / and B in favour of the decay constant measured on the lattice and the 
masses of the particles on the lattice; this removes any lattice spacing dependence arising from 
double trace operators. We conclude that there can never be any lattice-spacing dependence arising 
from NLO operators in mesonic scattering amplitudes, expressed in lattice-physical parameters, at 
next-to-leading order. 

Let us consider some simple examples to help make these arguments more concrete. A list 
of the NLO operators appearing in the staggered chiral Lagrangian is given in Each of these 
operators also appears in the NLO Lagrangian describing mixed action chiral perturbation theory 
with staggered sea quarks, with the rule that the taste matrices £ of Ref. [^] are only non-zero on the 
sea quark subspace; this is equivalent to replacing £ — > P s TjP s , where P s is the projector to the sea 
subspace, in spurions associated with taste violation. For example, one operator in the Lagrangian 
is 

ffi = ah\x [d^d^ 5 P s I?P£ 5 P s LP s ] . (2.3) 

We only want to use this operator at tree level, so we should take the sea quarks to vanish; this 
sets PgLPg — ► 0. Consequently, G\ does not contribute at NLO. For another example, consider the 
double trace operator 

G 2 = ahtx [d^d^L] str [% 5 P s tf P£ 5 P S LP S ] . (2.4) 

In this case, it helps to remember that £ is given in terms of the matrix of mesons <I> by £ = 
exp[2i<I>//], so when we set the sea quarks to zero the trace over the sea subspace is just a constant, 
namely the number of sea quarks N s . Thus, for the purposes of a next-to-leading order calculation 
@2 reduces to 

G 2 -> a 2 str [d^d^L] . (2.5) 

Notice that the effect of &i is to introduce an a 2 shift to the leading order parameter /. 

A similar argument applies to potential counterterms involving sea quark masses. These are 
associated with spurions arising from the sea sector and therefore they may not contribute to scat- 
tering amplitudes except to renomalise the parameters / and B; this dependence is removed upon 
switching to lattice-physical parameters. However, there is another, perhaps more physical, way 
of seeing why the sea quark masses do not contribute to the scattering. To understand this, we 
must digress briefly on %% scattering in SU (3) chiral perturbation theory. Evidently, the strange 
quark mass m s is a parameter of SU (3) #PT, and so one might think that the KK scattering length 
includes analytic terms involving m s . However, if we suppose that the strange quark is heavy and 
integrate it out, we must then recover SU (2) #PT. Now, chiral symmetry forces any m„ dependence 
in the analytic terms of the %% amplitude to occur in the form mr % rrv\. But the only counterterm in 



the on-shell SU (2) scattering amplitude (Eq. (2.2) with Apq = 0) is proportional to m\. It is not 



4 



GW Meson Scattering on a Staggered Sea 



Donal O'Connell 



possible to absorb m n m K into m n , so there can be no m s dependence in the SU (3) %% scattering 
amplitude. This is indeed the case [18]. 

Now, let us return to the situation with the valence-sea meson masses. For the purposes of 
this discussion, we can ignore the flavour-neutral and ghost sectors, reducing our theory from an 
SU(6\3) theory to an SU(6) theory. The sea quark dependence of this Si/ (6) chiral perturbation 
theory is analogous to the m s dependence of SU (3) #PT. A similar decoupling argument tells us 
that the sea quark masses cannot affect processes involving the valence sector provided one uses 
the analogues of on-shell parameters which are the lattice-physical parameters. We conclude that 
there can be no analytic dependence on the sea quark masses in a mesonic scattering amplitude. 
These rather abstract arguments have been verified by explicit computation Jl9[]. 

To conclude this section, we have shown so far that there can be no dependence on the lattice 
spacing a or the sea quark masses arising from the NLO operators in the MA^PT Lagrangian at 
next-to-leading order in the chiral power counting. However, the scattering amplitudes can still 
depend on these quantities through loop corrections from the LO Lagrangian. Note that since there 
is no counterterm available to absorb divergences proportional to the lattice spacing or the sea 
quark masses; this is a restriction on the form of the loop corrections. In favourable cases this 
restriction is strong enough to remove all dependence on the valence-sea mesons as occurred in the 
Tin amplitude; this removes all dependence on the unknown parameter Cmix- In the next section we 
will discuss individual mesonic processes and their structure. 



3. Meson meson scattering 

In a forthcoming work fll9| ] we present explicit expressions for the KK and K% scattering 
amplitudes at threshold. Here, we will restrict our comments to the dependence of these amplitudes 
on the valence-sea mesons. First, let us introduce some notation. It is useful to quantify the partial 
quenching in the mesonic sector with the parameters 

Aj u = rhjj — = 2B(ntj — m u ) + a 2 Ai (3.1) 
A rs = m rr — m ss = 2B(m r — m s ) + a 2 A/. (3.2) 

It is common to consider quantum theories which have unitary low energy sectors but which violate 
unitarity above some cutoff A. For example, 4 theory with a Pauli-Villars regulator is simply 
a theory or a light self-interacting scalar with a heavy ghost field which violates unitarity. The 
parameters A which we have introduced are measures of unitarity violation in the low-energy sector 
of QCD; when these parameters are zero, for example, the effective theory in Minkowski space has 
cuts at the points expected by the optical theorem. From the point of view of the low energy theory, 
one would like to tune these parameters to zero to remove as much of the unitarity violation as 
possible. 2 Note that this involves tuning the sea quarks to be lighter than the valence quarks; the 
high energy theory of quarks is non unitary but this is only seen above the non-perturbative scale 
of QCD, when perturbative computations are possible. With these ideas in hand, let us consider the 
mesonic scattering amplitudes. 

2 At higher orders, the A parameters will receive further corrections, but can generally be defined as the difference 
between lattice-physical sea-sea meson masses and lattice-physical valence-valence masses. 
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The KK scattering amplitude in lattice-physical parameters does not depend on the valence- 
sea mesons and hence the unknown parameter Cmix- When A„ = the expression for the scattering 
amplitude is 

_ 4m| 56m 4 K _ %m\ f m|\ / 10w|w| t 2m\ \ f m| 
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where the functions are known functions [19]. Note that the counterterm ^gl(m) in this expres- 



sion is the same as the counterterm in the physical scattering amplitude in on-shell parameters as 
we expect from the discussion above. 

The Kit scattering amplitude has a more complicated form than either the Tin or KK scattering 
amplitudes. This arises from the fact that the n and K meson masses are not equal, which induces 
a net momentum flow through loop propagators in the w-channel. Consequently, in this case the 
scattering amplitude does depend on the valence-sea mesons and therefore on Cmix- The part of the 
K% amplitude which depends on the sea quarks F = j, / , r is 

m K m % Y^ ( CVrflog — ™ ~ Cfs log — y - 2m K m 7t J(m 2 Fd ) +Am K m n (3.4) 



vs (4nf 2 ) 2 
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Although the amplitude does depend on the valence sea-mesons, it is straightforward to check that 
the \i dependence of the amplitude does not depend on the these mesons, in agreement with our 
comments above. Of course, the counterterm in the full K% amplitude in lattice-physical parameters 
is given by the QCD expression in on-shell parameters. 

4. Conclusions 

It is well-known that simulations with Ginsparg-Wilson valence quarks have no lattice spacing 
dependence at tree level of chiral perturbation theory as a consequence of the chiral properties of 
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GW quarks. We have shown that this persists at next-to-leading order in the chiral expansion 
for mesonic scattering processes, expressed in lattice-physical parameters, in the sense that the 
unknown analytic term which must be determined from the lattice, does not depend on the lattice 
spacing. In addition, we have shown that there is no unknown analytic terms involving the sea 
meson masses. This situation is favourable because it simplifies the extraction of the interesting 
physical information from a lattice simulation of these processes: in principle, one need only use 
one lattice spacing, for example. 

In addition, we have commented on the dependence of the KK and Kn scattering amplitudes 
on the unknown parameter Cmix which appears in the leading order chiral Lagrangian describing 
mixed action processes. The KK scattering amplitude does not depend on this parameter, but 
the Kn amplitude does; therefore, correct chiral extrapolation of Kn scattering data requires the 
determination of this parameter. 
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